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Abstract. We take a fresh look at the abundance of nitrogen in damped Lyman a systems (DLAs) with oxygen abundances 
between ~1/10 and ~ 1/100 of solar. This is a metallicity regime poorly sampled in the local universe and where QSO absorbers 
may hold clues to both the nucleosynthetic origin of N and the chemical evolution of high redshift galaxies. We combine new 
VLT UVES observations with others from the literature to form a sample of 10 DLAs in which the abundances of N and 
of one of the two a-capture elements O or S have been measured. The sample consists exclusively of high resolution, high 
signal-to-noise ratio data obtained with 8-10 m class telescopes. We confirm earlier reports that the (N/O) ratio exhibits a larger 
range of values than other ratios of heavy elements in DLAs; however, all 10 DLAs fall in the region of the (N/O) vs. (0/H) 
plot delimited by the primary and secondary levels of nitrogen production. Our results provide empirical evidence in support 
of the proposal that intermediate mass stars (4 ;S M/Mq ^ 7) are the main source of primary nitrogen, so that its release 
into the interstellar medium lags behind that of oxygen, which is produced by Type II supemovae. A high proportion (40%) 
of the DLAs in our sample have apparently not yet attained the full primary level of N enrichment; this finding may be an 
indication that the metallicity regime we are studying preferentially picks out galaxies which have only recently condensed out 
of the intergalactic medium and begun forming stars. Alternatively, the delay in the release of N following an episode of star 
formation may increase with decreasing metallicity, if stars of progressively lower masses than 4 Mq can synthesize N in their 
hotter interiors. In this general picture, the uniform value of (N/O) measured locally in metal-poor star-forming galaxies implies 
that they are not young, and is consistent with the presence of older stellar populations revealed by imaging studies with the 
Hubble Space Telescope. 

Key words. Cosmology: observations — Galaxies: abundances — Galaxies: ISM — Quasars: absorption lines 



1. Introduction 

The nucleosynthetic origin of nitrogen continues to be a sub- 
ject of considerable interest and discussion. There is general 
agreement that the main pathway is a six step process in the 
CN branch of the CNO cycle which takes place in the stel- 
lar H burning layer, with the net result that ^''N is synthesised 
from ^^C and ^^O. The continuing debate, however, centres 
on which range of stellar masses is responsible for the bulk 
of the nitrogen production. A comprehensive reappraisal of the 
problem was presented by Henry, Edmunds, & Koppen (2000) 
who compiled an extensive set of abundance measurements and 
computed chemical evolution models using published yields. 
Briefly, nitrogen has both a primary and a secondary compo- 
nent, depending on whether the seed carbon and oxygen are 
those manufactured by the star during helium burning, or were 
already present when the star first condensed out of the inter- 
stellar medium (ISM). 



/. /. Nitrogen and Oxygen in l-l II regions 

Observational evidence for this dual nature of nitrogen is pro- 
vided mainly from measurements of the N and O abundances in 
H II regions. (For simplicity in this paper we use parentheses to 
indicate logarithmic ratios of number densities; adopting the re- 
cent reappraisal of solar photospheric abundances by Holweger 
(2001), we have (N/U)q = -4.07; (O/H)© = -3.26; and 
(N/O)0 = -0.81). In H II regions of nearby galaxies, (N/O) 
exhibits a strong dependence on (O/H) when the latter is greater 
than ~ 2/5 solar; this is generally interpreted as the regime 
where secondary N becomes dominant.' At low metallicities 
on the other hand, when (O/H) < -4.0 (that is, < 1/5 solar), 



' As an aside, Henry et al. (2000) pointed out that the rise in (N/O) 
with (O/H) is steeper than would be normally expected for a purely 
secondary element and proposed that the effect is augmented by a de- 
creasing O yield with increasing metallicity (Maeder 1992). 
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Table 1. Damped Lya Systems Observed 



QSO 


V 






N{K 1)" 


Exp. Time 


Resolution 


Wavelength Range 






(mag) 






(cm-2) 


(s) 


(kms-^) 


(A) 




Q1409+095 


18.6 


2.856 


2.45620 


3.5 X 10^° 


16200 


7 


3793^989, 6726-10255= 


15-25 


Q1444+014 


18.5 


2.206 


2.08681 


1.6 X 10^° 


18 000 


7 


3295-6650"* 


10-20 


Q2206-199 


17.33 


2.559 


2.07623 


2.7 X 10^° 


24300 


7 


3295-10 255'' 


20-25 


' Determined by profile fitting to the damping wings of the Lyman a absorption line (see section 2 and Figure 1). Typical error 



mN(H I) is ±10%. 

^ The signal-to-noise ratio varies along each spectrum. These values (per pixel) refer to the QSO continuum in the regions of 

the N I AA1200 and O I A1302.2 lines. 
° With a wavelength gap between 8520 and 8667 A. 
With some wavelength gaps. 



N is mostly primary and tracks O; this results in a plateau at 
(N/O) ~ -1.5. 

The principal sources of primary N are thought to be in- 
termediate mass stars (4 £ M/M© ^ 7) during the asymp- 
totic giant branch (AGB) phase. Henry et al. (2000) showed 
that, integrating along the IMF the N yields from intermediate 
mass stars by van den Hoek & Groenewegen (1997) and the 
O yields from massive stars by Maeder (1992) at metaUicity 
Z = 1/20.^0, one obtains (N/O) = —1.41 in good agreement 
with the observed plateau. 

A corollary of the hypothesis that intermediate mass stars 
are the main producers of primary nitrogen is that its release in 
the ISM should lag behind that of O, since the latter is widely 
believed to be produced by massive stars which explode as 
Type II supernovae (SN) soon after an episode of star forma- 
tion. Henry et al. (2000) calculated this time delay to be approx- 
imately 250 Myr; at low metaUicities the (N/O) ratio could then 
perhaps be used as a clock with which to measure the past rate 
of star formation, as proposed by Edmunds & Pagel (1978). 
Specifically, in metal-poor galaxies which have only recently 
experienced a burst of star formation one may expect to find 
values of (N/O) below the primary plateau at (N/O) ~ —1.5. 

This scenario has been discussed at length over the last few 
years with claims of both conflicting and supporting evidence 
from the data. Izotov & Thuan (1999) and Izotov et al. (2001) 
were struck by the constant values of all the element ratios, 
including (N/O), which they measured in blue compact dwarf 
(BCD) galaxies with (O/H) ^ — 4.4 (corresponding to an oxy- 
gen abundance less than ~1/15 of solar). This finding led them 
to the radical proposal that (a) these are galaxies undergoing 
their first burst of star formation, and (b) all the elements they 
observed have a primary origin in massive stars, so that they are 
released at the same time as O, thereby disposing altogether of 
the notion of a time delay. On the other hand, a variety of stud- 
ies using high resolution Hubble Space Telescope (HST) im- 
ages have shown that many of these BCDs host old stellar pop- 
ulations, with ages greater than 1 Gyr (e.g. Schulte-Ladbeck et 
al. 2001 and references therein; Crone et al. 2002). Generally, 
these low mass galaxies tend to have low rates of star formation 
(even though they are observed during a bursting phase). If the 



lag between O and N production is only 250 Myr, as proposed 
by Henry et al. (2000), it is perhaps not surprising that most 
of them lie near the primary value of the (N/O) ratio expected 
from intermediate mass stars (Pilyugin 1999). 

In contrast, a recent survey by Contini et al. (2002) found 
that UV-selected galaxies at intermediate redshifts (0<z<0.4) 
exhibit a wide range of (N/O) values at a given (O/H), in many 
cases well below the primary level. While Contini et al. inter- 
pret their results as evidence in favour of a delayed produc- 
tion of primary N, one may well wonder at the high fraction of 
galaxies apparently caught within 250 Myr since the last major 
episode of oxygen enrichment. Possibly the UV-continuum se- 
lection technique strongly favours galaxies with high rates of 
star formation in the recent past. 

1.2. Nitrogen and Oxygen in QSO Absorption Line 
Systems 

As pointed out by Pettini, Lipman, & Hunstead (1995), clues 
to the nucleosynthetic origin of nitrogen can also be provided 
by QSO absorption line systems, particularly the high column 
density (iV(H I) > 2 x 10^° cm~^) damped Lyman a absorbers 
(DLAs). These are thought to represent an early stage in the 
evolution of galaxies, when most of their baryonic mass was 
in the interstellar medium. Apart from the obvious interest in 
measuring element abundances in the distant past, when galax- 
ies were young, one of the advantages of DLAs is that they 
are generally of low metaUicity, approximately between 1/10 
and 1 /lOO solar (Pettini et al. 1999; Prochaska & Wolfe 2002). 
Thus, they probe a regime where local H II region abundance 
measurements are sparse or non-existent and where the effect 
of a delayed production of primary nitrogen should be most 
pronounced. 

There are practical difficulties, however, in measuring 
(N/O) in DLAs. The resonance Unes most easily accessible, 
the N I triplet AA1199.5, 1200.2, 1200.7and O I A1302.2, have 
very different optical depths. While the latter is normally sat- 
urated in DLAs, the former is weak and can be difficult to de- 
tect, particularly as it falls within the Lyman a forest. Pettini 
et al. (1995), and subsequently Lu, Sargent, & Barlow (1998) 
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and Centurion et al. (1998), attempted to circumvent the first 
problem by using the abundances of other a-capture elements, 
such as Si and S, as proxies for (O/H). These studies found that 
the (N/a) ratio in DLAs spans a wide range of values, broadly 
within the range bracketed by the predictions for primary and 
secondary production of nitrogen. However, these results have 
been criticized by Izotov & Thuan (1999) and more recently 
Izotov, Schaerer, & Charbonnel (2001). Unlike O and N, which 
are mostly neutral in H I gas. Si and S are singly ionised. 
Since in principle Si II and S II can also occur in H II regions, 
Izotov and collaborators have proposed that unaccounted ioni- 
sation corrections, rather than real abundance variations, are the 
cause of the observed scatter of the (N/a) ratios in DLAs. Their 
model, however, makes the radical assumption that essentially 
all the metals are located in (presumably self-enriched) H II re- 
gions, and that the metalUcity of the neutral phase is negUgible. 

1.3. Motivations for tfie Present Study 

In an attempt to clarify this confusing state of affairs, we have 
begun a new programme of observations with the Ultraviolet- 
Visual Echelle Spectrograph (UVES) on the VLT (Kueyen) 
telescope (Dekker et al. 2000). In order to minimize the prac- 
tical difficulties outlined above, we targeted DLAs which 
(a) have relatively low values of hydrogen column density 
(iV(H I) ^ 4 X 10^" cm~^), thus increasing the probabihty that 
the O I A1302.2 line may not be strongly saturated), and (b) 
are at relatively low redshifts (zabs ^ 2.7), where the Lyman a 
forest begins to thin out. Despite this observing strategy, which 
capitalises on the superior performance of UVES at ultraviolet 
wavelengths (D'Odorico et al. 2000), our endeavours have been 
only partially successful, as we shall see. In addition, we bring 
together a number of relevant abundance measurements which 
have been published since the surveys by Lu et al. (1998) and 
Centurion et al. (1998) for a comprehensive reappraisal of the 
abundances of oxygen and nitrogen in DLAs. 

2. Observations and Data Reduction 

Details of the DLAs observed are collected in Table 1. The 
data were secured over the three nights of 28-30 May 2000. 
By using different cross-disperser gratings and dichroic filters 
we recorded the spectra of three QSOs over the wavelength 
ranges hsted in the penultimate column of Table 1. With a 1 
arcsec wide entrance slit and 2x2 binning on the CCDs, the 
resulting spectral resolution was between 6.7 and 7.1km s~^ 
FWHM, measured from the widths of emission Unes of the Th- 
Ar hollow-cathode lamp used for wavelength calibration. The 
observations consisted of a number of one hour long exposures; 
the QSOs were moved along the slit between exposures so as 
to use different portions of the detectors to record the spec- 
trum of each object. The sUt was aUgned at the parallactic angle 
throughout the observations. 

The two-dimensional images were processed with the 
UVES standard pipeline software; see for example Ellison, 
Ryan, & Prochaska (2001) for a description of the main steps. 
The individual spectra extracted for each QSO were mapped 
onto a common, hnear, vacuum heUocentric, wavelength scale; 



the interval between successive wavelength bins was set to 
the average value of the original CCD pixels, typically i5A = 
0.03 A and 0.04 A in the blue and red wavelength regions re- 
spectively. The individual extractions were then co-added, with 
weights proportional to the S/N ratios, to produce final blue and 
red spectra for each QSO. With several independent exposures 
at our disposal, we were able to exclude from the sum wave- 
length bins which deviated from the mean by several standard 
deviations, because of residuals from the subtraction of cosmic- 
ray events, strong sky emission lines, and other invalid pixels. 

Figure 1 shows portions of the final spectra in the re- 
gions near the damped Lyman a lines at ^abs = 2.45620 in 
Q1409+095, Zabs = 2.08681 in Q1444+014, and Zabs = 
2.07623 in Q2206-199, after division by the underlying QSO 
continuum. Superposed on the data are theoretical damped ab- 
sorption profiles corresponding to the values of N(li I) listed 
in Table 1 ; the typical error in the determination of iV(H I) is 
about 10%. Our observations revealed the existence of a sec- 
ond high column density absorption system in Q1409+095, 
at ^abs = 2.66820. For completeness, we have included it in 
Figure 1; the profile fit shown is for A/'(H I) = 5 x 10^^ cm'^. 

3. Metal Absorption Lines in tlie Four DLAs 

3.1. Rationale 

As is the case with modem echelle spectra, our observations of 
the three QSOs in Table 1 provide a wealth of information on 
several absorption systems, as well as on the Lyman a forest. 
In the present study we are primarily interested in the abun- 
dances of nitrogen and oxygen, detected respectively via the 
N I AA1199.5, 1200.2, 1200.7 triplet and the O I A1302.2 line. 
Also of relevance to our analysis are the abundances of siU- 
con (measured via Si II A1304.4) and iron (several absorption 
lines of Fe II throughout the spectral regions covered). Silicon 
is an a-capture element which can provide an additional check 
on the oxygen abundance in cases where the O I A1302 line is 
saturated and the corresponding value of N{0 I) is uncertain. 
Sulphur would also be useful in this context, but unfortunately 
the relevant absorption lines are either blended or undetected 
in most of the cases considered here. Iron is of interest because 
it is thought to be produced primarily by Type la supemovae; 
like nitrogen, its release into the interstellar medium should be 
delayed relative to that of oxygen. The four elements in ques- 
tion, N, O, Si, and Fe, respond differently to depletion onto 
dust grains. In the local ISM, N and O are normally present 
in the gas phase with near-solar abundances; Si and Fe, on the 
other hand, are susceptible to dust depletion, the latter normally 
more than the former (e.g. Savage & Sembach 1996). Thus, by 
considering all four elements together, it may be possible to 
decouple nucleosynthesis and depletion effects. 

3.2. Profile Fitting and Abundance Determinations 

Normalised profiles of the absorption lines mentioned above 
were fitted with Voigt profiles using the VPFIT package.^ 

^ VPFIT is available at http://www.ast.cam.ac.uk/~rfc/vpfit.html 
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Table 2. Component Structure of the Metal Absorption Lines 





Component 
Number 


■^abs 


L 



(kms-^) 


log iV(U Ij 

(cm-2) 


log 11) 
(cm-2) 


log i V (re 11 ) 

(cm-2) 


log iV(JN 1) 

(cm-2) 


Q1409+095 


1 


2.45597 


12.9 


14.76 ±0.1 


13.82 ± 0.03 


13.46 ± 0.03 


< 12.80 




2 


2.45630 


5.8 


14.27 ±0.03 


13.08 ±0.03 


12.81 ± 0.03 


< 12.50 




3 


2.45648 


7.3 


14.82 ±0.1 


13.63 ±0.03 


13.29 ±0.03 


< 12.77 


Q 1409+095 


1 


2.66820 


7.5 


15.27 ±0.1 


14.02 ± 0.03 


13.85 ±0.03 


< 13.45 


Q2206-199 


1 


2.07623 


6.1 


15.24 ±0.09 


13.65 ± 0.02 


13.37 ± 0.02 


< 12.88 



VPFIT deconvolves the composite absorption profiles into the 
minimum number of discrete components and returns for each 
the most likely values of redshift z, Doppler width b (km s~^), 
and column density N (cm~^) by minimizing the difference 
between observed and computed profiles. The profile decom- 
position takes into account the instrumental point spread func- 
tion of UVES. As we shall see shortly, the velocity structure of 
three out of the four DLAs studied here is very simple. In each 
case the same set of values of z and b was found to provide a 
satisfactory fit to the absorption lines of O I, Si II, and Fe II. 
Oscillator strengths and rest wavelengths were adopted from 
the up to date compilation of atomic parameters maintained by 
J. X. Prochaska at http://kingpin.ucsd. edurhiresdla/atomic.dat 
(see also Prochaska et al. 2001). The solutions returned by 
VPFTT are collected in Table 2. 

The error estimates tT(log N) returned by VPFIT were typ- 
ically < 0.05 dex for the optically thin Si II and Fe II lines and 
< 0.1 dex for the saturated O I lines. To test whether these er- 
rors are realistic we performed Monte Carlo type simulations 
similar to those described by Bowen, Blades, & Pettini (1995) 
and Bowen, Pettini, & Blades (2002). Briefly, for each DLA 
we used as the starting point the theoretical spectrum produced 
by VPFIT with the best fit parameters given in Table 2. From 
this template we generated 300 spectra with the same S/N as 
the original data; these were then refitted to deduce values of 
b and A''. From gaussian fits to the distributions of 300 values 
of b and N we determined values of aib) and cr(Iog N). The 
former are typically ±0.2 km s^^; the latter, which are com- 
parable to the values returned by VPFIT, are listed in Table 2. 
Despite the good agreement between the two methods, it must 
be remembered that these errors reflect primarily the S/N of the 
data and do not include the possibility that the velocity disper- 
sion parameter b may be different between the different ions 
considered. The consequences of letting b vary from ion to ion 
are discussed below for each DLA individuaUy. 

The species observed are the major ionisation stages of 
their respective elements in H I regions. Thus, to obtain their 
abundances relative to hydrogen we simply divide the column 
densities returned by VPFIT by iV(H I) (after adding together 
the values for individual velocity components within the same 
DLA, when more than one component is present). Finally these 
values are referred to the solar abundance scale (Holweger 



2001). The results are collected in Table 3, where the solar ref- 
erence values are listed in the footnotes. 

This procedure would overestimate the abundances of Si 
and Fe if some of the Si II and Fe II absorption we detect were 
due to ionised gas associated with the neutral DLA. Conversely, 
the abundances of O and N would be underestimated if these 
elements were partly ionised in the H I region. In the cases un- 
der consideration, however, there is httle evidence that these 
reservations are justified. First, all four species exhibit very 
similar velocity structure; we see no indication of components 
present in the first ions and missing in the neutrals, or vice 
versa. Second, the recent reappraisal of this problem by Vladilo 
et al. (2001) shows that at the values of N(K I) of the DLAs 
studied here, the ionisation corrections to the abundances de- 
termined as above are smaller than the typical error of ±0.15 
dex of our abundance estimates. We now briefly comment on 
each DLA in turn. 

3.3. The 2abs = 2.45620 DLA in Q1409+095 

The Si II A1304.4 and Fe II A2344.2 lines in this DLA are un- 
saturated and are weU reproduced with three absorption com- 
ponents spanning a velocity range Av = 44km s~^ (see Figure 
2 and Table 2). This is a very metal-poor system. Adding to- 
gether the column densities for the three components we de- 
duce that the abundance of Si is 1/100 of solar, and that of Fe 
only little more than 1/200 of solar (see Table 3). 

Even at these low abundances, and with a hydrogen col- 
umn density of 3.5 x lO^^cm"^ (near the lower limit of 
2 X 10^° cm~^ commonly adopted as the definition of a DLA), 
the O I A1302.2 line is saturated (see Figure 2). If we adopt 
the same set of absorption parameters as for Si II and Fe II, we 
deduce an oxygen abundance of 1/135 of solar, in reasonable 
agreement witii that of Si. The good model fit to O I A1302.2 
makes it unlikely that we have overestimated N(0 I) (the line 
would be broader than observed). The oxygen abundance could 
be higher if our model underestimates the saturation of the 
O I A1302.2 line; this possibility, however, does not alter the 
interpretation of our results below because it would tend to em- 
phasize the 'anomalously' low (N/O) ratio (see §5). 

The N I A1200 triplet at Zabs = 2.45620 falls within the 
Lyman a forest of this Zgm = 2.856 QSO. Figure 2 shows 
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QSO 




log 7V(H 1)" 
(cm-2) 


log iV(0 1) 
(cm-') 


log iV(Si II) 


log A'^CFe II) 
(cm-2) 


log Af(N I) 
{cm-') 


[OlYLf 


[Si/H]'' 




[N/H]*' 


Q1409+095 


2.45620 


20.54 


15.15 ±0.06 


14.08 ± 0.02 


13.74 ± 0.02 


< 13.19 


-2.13 


-2.00 


-2.25 


< -3.28 


Q1409+095 


2.66820 


19.70 


15.27 ±0.1 


14.02 ± 0.03 


13.85 ± 0.03 


< 13.45 


-1.17 


-1.22 


-1.30 


< -2.18 


Q2206-199 


2.07623 


20.44 


15.24 ±0.09 


13.65 ± 0.02 


13.37 ± 0.02 


< 12.88 


-1.94 


-2.33 


-2.52 


< -3.49 



Typical error in log Niii I) is ±0.04. 

In the usual notation: [X/H] = log(X/H) — log(X/H)0. Solar abundances are from the recent updates by Holweger (2001): 
log(O/H)0 = -3.264; log(Si/H)0 = -4.464; log(Fe/H)0 = -4.552; andlog(N/H)Q = -4.069. 



the best fit to this spectral region obtained with the model 
parameters of the other three lines, but it is difficult to be 
confident that the weak absorption we see, particularly in the 
A1200.2 transition, is actually N I and not a Lyman a for- 
est line. Consequently, we prefer to adopt the corresponding 
7V(N I) as an upper Umit to the column density of N I, from 
which we deduce that the abundance of N is less than 1/2000 
of solar. 

3.4. The Zabs = 2.66820 sub-DLA in Q1 409+095 

There is a second high column density system in Q1409-I-095; 
although its N(}i 1)— 5.0 x 10^^ cm~^ is below the Umit nor- 
mally adopted for DLAs (such systems are sometimes referred 
to as sub-DLAs), the damping wings of the Lyman a Une are 
easily recognisable in the UVES spectrum (see Figure 1) and 
7V(H I) is well determined. 

The O I, Si II, and Fe II absorption Unes in this system 
are well reproduced with a single absorption component with 
a narrow velocity dispersion, b = 7.5 km s~^ (see Figure 3); 
although there may be further substructure within this com- 
ponent, it is difficult to resolve it unambiguously in our data. 
All three lines give consistent abundances of between 1/15 and 
1/20 of solar (see Table 3). Again, the N I triplet is blended 
with Lyman a forest absorption. While there are features which 
match the expected positions of two of the N I lines, we take 
the fit shown in Figure 3 as indicative of the maximum amount 
of N I which could be present. In this case, N is less abundant 
than 1/150 of solar. 

3.5. The zabs = 2.08681 DLA in Q1444+014 

Unlike the other DLAs in this study, this system exhibits a very 
complex velocity structure. The fit to the Si II A1304.4 line 
requires at least 13 components spread over the velocity in- 
terval Av ~ 350 km s~^, from -170 to +176km s~^ rela- 
tive to Zabs ~ 2.08681 (see Figure 4). Unfortunately, many of 
the components are strongly saturated in O I; while the weaker 
ones are well matched by the model parameters returned by 
VPFTT for Si II, we cannot determine a reliable value of the 
total column density of O I as it is dominated by the heavily 



saturated absorption between w ~ and —150km s~ . In the 
bottom panel we show the partial fit, by the weaker components 
only, to the region encompassing the N I triplet. The combina- 
tion of 13 absorption components which partially overlap in 
the closely separated triplet fines, together with blending with 
Lyman a forest lines, makes it impossible to estimate the abun- 
dance of N. As we are unable to determine (N/O) in this system, 
we do not consider it further in this paper. 

3.6. The Zabs = 2.07623 DLA in Q2206- 199 

This is a well-known DLA, one of the first to be studied at 
echelle spectral resolution thanks to the brightness of the QSO 
(Pettini et al. 1990; Pettini & Hunstead 1990). It has a re- 
markably simple velocity structure — the absorption apparently 
arises in a single component with 6 ~ 6 km s^^ (Figure 5 and 
Table 2). This property makes it one of the few DLAs suitable 
for measuring the abundance of deuterium (Pettini & Bowen 
2001). It also exhibits a very low metallicity; O, Si and Fe are 
below solar by factors of 90, 210, and 330 respectively (Table 
3). For the elements in common (Si and Fe), our results are in 
good agreement (within 0. 1 dex) with the Keck data analysed 
by Prochaska & Wolfe (1997). As can be seen from Figure 5, 
we have a tentative detection of N I from which we deduce an 
upper limit to the abundance of N of only 3 x 10"'' of solar. 

It is interesting that in this system we deduce [0/SiJ = 
+0.4: ± 0.1 (Table 3). Possibly, this is a hint of the trend of 
increasing [O/Si] with decreasing metallicity seen in Galactic 
halo stars with [Fe/H| ^ — 2 (Edvardsson et al. 1993; Ryan, 
Norris, & Beers 1996; Israelian et al. 2001). Less interesting 
explanations are that Si is depleted onto dust (although a deple- 
tion by a factor of 0.4 in the log seems somewhat high, given 
the low metallicity of this DLA — see Pettini et al. 1997) or 
that we have overestimated N{0 I) when fitting the saturated 
A1302 absorption line. The abundance of S, an undepleted ele- 
ment which tracks O more closely than Si (Israelian & Rebolo 
2001; Takada-Hidai et al. 2001 — see discussion in §4) would 
help discriminate between these three possibihties, but unfor- 
tunately all three transitions of the S IIA1256 multiplet in this 
system are blended with Lyman a forest lines. 
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4. The DLA Sample 

In the last few years there have been several new measurements 
of the abundance of N, generally reported as part of abundance 
studies of individual DLAs. We have collected these data in 
Table 4 where references to the original papers can also be 
found. All of the measurements in the table are from spec- 
tra of high signal-to-noise ratio and high (echelle) resolution 
obtained with 8- 10 m class telescopes. For simplicity we hmit 
ourselves to absorption systems with NiH I)> 2 x 10^° cm~^; 
inclusion of similar data available for a few sub-DLAs would 
not change any of our conclusions. 

The full sample consists of 10 DLAs in which the abun- 
dance of N has been measured (seven detections and three 
upper limits). In one half of the cases the abundance of O is 
available directly, either because the O I A 1302. 2 line is not 
strongly saturated, or from weaker, unsaturated, transitions to 
higher energy levels. In the other five cases we take N{S II) as a 
proxy for A'^(0 1), after correcting for the solar abundance ratio 
(O/S)0 = -1-1.54 (Grevesse & Sauval 1998; Holweger 2001). 
The justification for this approach is as follows. 

First, in Galactic stars the (O/S) ratio appears to be ap- 
proximately solar over a wide range of metallicities down to 
[O/H] < - 2 (IsraeHan et al. 2001; Israehan & Rebolo 2001; 
Takada-Hidai et al. 2001; Nissen et al. in preparation). Second, 
in the diffuse Galactic interstellar medium both elements show 
little affinity for dust (Savage & Sembach 1996); dust correc- 
tions for O and S are almost certainly unimportant for our pur- 
poses at the reduced levels of dust depletions typical of metal- 
poor DLAs (Pettini et al. 1997; Vladilo 1998). Finally, differ- 
ential ionisation corrections between O I and S II are likely 
to be small. Theoretically, the photoionisation calculations by 
Vladilo et al. (2001) show that, for all the models they consid- 
ered, if log (O/S) = log [iV(0 l)/N(S II)] -I- C, then |C| < 0.2 
( < 0.1) when log A''(H I)> 20.2 (> 20.6). This condition 
is met by 5 (4) out of the 5 DLAs in Table 4 for which we 
use S as a proxy for O. For all the DLAs considered here 
|C|[N/0] < 0.002 so that log (N/O) =log [NiN l)/NiO I)]. 

Some authors (e.g. Lu et al. 1998) have also used Si 
to deduce the abundance of O, by assuming that in DLAs 
A'^(0 T)/N(Si II) = (O/Si)0. This assumption, however, is less 
secure than the equivalent one for S and has been criticised, 
for example, by Matteucci, Molaro, & Vladilo (1997). First, 
Si can be depleted onto dust. Second, it is unclear whether the 
abundances of O and Si do track each other at low metallicities 
(Edvardsson et al. 1993; Ryan, Norris, & Beers 1996; Israelian 
et al. 2001); theoretically they are not expected to, if some 
of the Si is produced by type la supemovae (Matteucci et al. 
1997). Neither of these complications appUes to S. In the past, 
the use of Si as an indicator of the O abundance had been mo- 
tivated by the paucity of measurements of [O/H] and [S/H] in 
DLAs; with the larger body of data now available this approach 
is no longer necessary. For this reason, we restrict ourselves to 
O and S in the following analysis. 



4. 1. Nitrogen and Oxygen in DLAs 

We are now ready to compare the abundances of N and O in 
high redshift DLAs with those measured in H II regions in the 
local universe (Figure 6). For the latter, we use the compilation 
assembled by Henry et al. (2000 — references to the original 
surveys are given in the figure caption), augmented by some 
additional work (Ferguson, Gallagher, & Wyse 1998; van Zee 
2000). The dashed lines in the figure are approximate bound- 
aries of the region in the (N/O) vs. (O/H) plot where we may 
expect DLAs to fall. The Une labelled 'Secondary' is simply an 
extrapolation to low (O/H) values of the (N/O) vs. (O/H) trend 
at high metallicities, where secondary N presumably domi- 
nates. If the extrapolation to low metalUcities is valid (we do 
not know that it is), this is the minimum amount of N (relative 
to O) which one may expect to find, irrespectively of the pre- 
vious star formation history of the system under consideration. 
The 'Primary' level is just the plateau at (N/0)= —1.5 indi- 
cated by existing measurements in low metallicity H II regions. 
We stress that these are just approximate empirical Umits which 
are not derived from a full chemical evolution treatment of the 
available data (e.g. Larsen, Sommer-Larsen, & Pagel 2001) 

As can be seen from Figure 6, the DLAs in our sample do 
indeed populate the region of the diagram delimited by the bro- 
ken lines. Six of the DLAs considered exhibit levels of N en- 
richment which are similar to those of the most metal-poor H II 
regions. On the other hand, in four cases N appears to be below 
the 'Primary' fine, by up to nearly one order of magnitude. 



5. Discussion 

The results presented here, in Figure 6 and Table 4, confirm 
the earher conclusion by Lu et al. (1998) that the (N/O) ra- 
tio in DLAs exhibits a larger range of values than that of any 
other pair of heavy elements which are not depleted onto dust 
(Prochaska & Wolfe 2002). We now consider possible causes 
of this scatter. 

The simplest explanation may be just measurement error. 
Although the errors on N{N I) and A''(0 I) quoted by the 
sources referenced in Table 4 are normally less than 0.1 dex 
(and always less than 0.2 dex), it must be remembered that the 
N I Unes (a) are generally weaker than most other metal transi- 
tions used for abundance work in DLAs, and (b) occur within 
the Lyman a forest where blending can complicate their analy- 
sis. It is thus conceivable that the true uncertainties, particularly 
in the N abundance, may have been underestimated. On the 
other hand, if measurement error were the sole cause of the dis- 
persion of the DLA points in Figure 6, it would be a remarkable 
coincidence that all 10 values in our sample should fall within 
the boundaries delimited by the Primary and Secondary Unes in 
the figure. In principle one would expect measurement errors to 
scatter points above the Primary line and, at the larger values 
of (O/H), below the Secondary boundary. For many of the low 
points in Figure 6, it would be very hard to understand how 
an absorption Une ~ 3 — 5 times stronger than observed could 
have been missed. We conclude, therefore, that measurement 
errors are not the principle cause of the dispersion observed. 



Max Pettini et al.: Nitrogen and Oxygen in Damped Lyman a Systems 7 



Table 4. N, O, S, and Fe Abundance Measurements in DLAs 



QSO 


^abs 


log N{n I) 


log iV(N 1) 


log N{0 I) 


log Ar(S II) 


(0/H)+12'' 


(N/O)" 


(Fe/O)" 


Ref.'' 






(cm-2) 


(cm-2) 


(cm-2) 


(cm-2) 










QOOOO-2620 


3.3901 


21.41 


14.73 


16.45 




7.04 


-1.72 


-1.58 


1 


QO 100+ 1300 


2.3090 


21.32 


15.29 




15.12 


7.34 


-1.37 


-1.56 


2, 3 


Q0201 + 1120 


3.38639 


21.26 


15.33 




15.21 


7.49 


-1.42 


-1.40 


4 


J0307-4945 


4.46658 


20.67 


13.57 


15.91 




7.24 


-2.34 


-1.70 


5 


Q0347-3819 


3.02486 


20.63 


14.89 


16.64 




8.01 


-1.75 


-2.21 


6 


Q1409+095 


2.45620 


20.54 


< 13.19 


15.15 




6.61 


< -1.96 


-1.41 


7 


GB 1759+7539 


2.62528 


20.76 


14.99 




15.21= 


7.99 


-1.76" 


-1.81 


8,9 


Q2206-199 


2.07623 


20.44 


< 12.88 


15.24 




6.80 


< -2.36 


-1.87 


7 


Q2343+1230 


2.4313 


20.34 


14.67 




14.71 


7.91 


-1.58 




2 


Q2348-1444 


2.27936 


20.56 


< 13.22 




13.73 


6.71 


< -2.05 


-1.48 


3 



^ When the oxygen abundance is not available, S has been used as a proxy for O by assuming the solar ratio (O/S)© = +1.54 

(Grevesse & Sauval 1998; Holweger 2001); see text (§4) for justification. 
^ References— 1: Molaro et al. (2000); 2: Lu, Sargent, & Barlow (1998); 3: Prochaska & Wolfe (1999); 4: Ellison et al. (2001); 

5: Dessauges-Zavadsky et al. (2001); 6: Levshakov et al. (2002); 7: This work; 8: Outram, Chaffee, & Carswell (1999); 9: 

Prochaska et al. (2002). 

° Prochaska et al. (2002) found that only about half of the metals in this DLA arise in a neutral component, while the reminder 
are in partially ionised gas. It is unclear what fraction of S in the ionised gas is S II; in the extreme case that half of the observed 
S n is in an ionised region where O I and N I are absent, the value of (N/0) for this DLA should be increased by +0.3 dex. 



We do not subscribe to the view proposed by Izotov et al. 
(2001) that unrecognised ionisation corrections are responsible 
for the low (N/O) values. As explained above, in about half 
of the cases (N/0)dla is comparable to the value measured 
in nearby metal-poor H II regions; there is no evidence that 
these DLAs are different from the rest of the sample in their 
level of ionisation. Similarly, there is no separation in Figure 

6 between DLAs where (O/H) has been measured directly and 
those where it has been deduced by reference to S. 

Some may be concerned about possible systematic differ- 
ences between abundances measured via absorption lines in 
cool interstellar H 1 clouds and those deduced from the analysis 
of nebular emission lines from H II gas. However, there are no 
precedents for such systematic differences in the solar vicinity 
(Meyer, Cardelli, & Sofia 1997; Meyer, Jura, & Cardelli 1998) 
and, in any case, here we are dealing with a scatter rather than 
a systematic offset. 

Qualitatively, our findings are consistent with, and indeed 
provide empirical evidence in favour of, an origin of primary 
nitrogen in intermediate mass stars. The observed range in the 
(N/O) ratio at low (O/H) is, in this picture, a natural conse- 
quence of the delayed release of N into the ISM relative to the 
O produced by Type II supemovae. However, there are puzzling 
aspects in this scenario too, in that the proportion of DLAs with 
values of (N/O) below the primary level seems at first sight 
to be surprisingly high. Recall that Henry et al. (2000) esti- 
mated the time delay between N and O release into the ISM to 
be ~ 250 Myr, if the main source are stars in the mass range 

7 — 4 Mq . The median redshift of the DLA sample in Table 
4 is (z) = 2.5; in a VIm = 0.3, = 0.7, h = 0.65 cosmol- 
ogy this corresponds to a look-back time of 11 . 7 Gyr. ^ There is 

^ The age of the universe in this cosmology is 14.5 Gyr. 



now evidence that star formation in the universe was already at 
high levels well before this epoch, certainly by 2: ^ 4.5 (Steidel 
et al. 1999), and possibly up to z ~ 6 (Shapley et al. 2001; 
Becker et al. 2001; Dawson et al. 2001). In this cosmology, 
z = & corresponds to a look-back time of ~ 13.5 Gyr. Thus, 
the time interval between z = 6 and our median (zdla) =2.5 
is seven times higher than the estimated time delay for the re- 
lease of primary nitrogen. If the galaxies giving rise to DLAs 
formed continuously between z = 6 and 2.5, naively one may 
expect to find only one DLA out of seven in the interim period 
when primary nitrogen has not yet caught up with the oxygen 
released by the same episode of star formation. Instead we find 
at least four such cases out of ten (see Figure 6). 

Perhaps we are just seeing the effects of small number 
statistics and a larger sample of (N/O) measurements will in 
future reveal that the present 40% is an overestimate of the 
number of DLAs with a N deficiency. However, if such a high 
fraction is confirmed by future observations, we can think of 
two explanations for this apparent puzzle. The four DLAs with 
values of (N/O) clearly below the primary level in Figure 6 
have abundances (O/H) less than ~ 1/30 of solar. Thus, they 
have presumably experienced little star formation up to the 
time when we observe them; in a closed box model, they would 
have turned less than 1/30 of their gas into stars. It therefore 
seems at least plausible that these DLAs arise preferentially in 
young galaxies, which have only recently condensed out of the 
intergalactic medium. In other words, the very fact that the dif- 
ference between primary and secondary nitrogen is most pro- 
nounced at low metallicities introduces a bias in our sample 
such that it may include a relatively high proportion of galax- 
ies caught in the initial stages of their chemical evolution (e.g. 
Cen et al. 2002). 
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A second possibility is that the time delay for the release 
of primary nitrogen may be longer when metal abundances 
are lower. This dependence may arise if at lower metallici- 
ties stars of progressively lower masses can synthesize and re- 
lease nitrogen, because their interiors are hotter (Lattanzio & 
Forestini 1999) and/or through the increasing importance of 
rotation (Meynet & Maeder 2002a,b; see also Marigo 2001). 
The slope of the initial mass function would then shift the peak 
of the N production to lower mass stars, with longer evolution- 
ary times. While these effects remain to be fully quantified, 
there have been suggestions that at the metallicities of DLAs 
the full release of primary N may lag behind that of O by up 
to 800 Myr (Lattanzio et al., in preparation). If confirmed, this 
proposal would go a long way towards explaining the high pro- 
portion of DLAs which are underabundant in N at (z) = 2.5. 

5. 1. Clues from the Abundance of Iron? 

It should be possible to test the hypothesis that the dispersion of 
(N/O) values in DLAs is a consequence of the delayed release 
of primary nitrogen by considering the abundances of other el- 
ements synthesized by low mass stars. Observationally, iron is 
one of the more easily accessible through several rest frame ul- 
traviolet transitions. Measurements of (Fe/H) are available for 
9 out of the 10 DLAs in our sample^; the corresponding values 
of (Fe/O) are hsted in Table 4. 

It is generally thought that approximately 2/3 of the iron is 
produced in Type la supernovae. While it is still unclear which 
types of binary system give rise to such events, there seems to 
be a consensus that the evolutionary timescales of Type la SN 
progenitors are of the same order as, or longer than, those of the 
7 — 4 Mq stars which are the source of primary nitrogen (e.g. 
Matteucci & Recchi 200 1 ) . From this it follows that DLAs defi- 
cient in N (in the sense of having less N than the primary level) 
should also exhibit sub-solar (Fe/O) ratios. We examine this 
point in Figure 7, by plotting the quantity [(N/O) + 1.5] (that 
is, the difference between the measured (N/O) and the primary 
plateau at (N/0)= -1.5) vs. [(Fe/O) + 1.29] (the difference 
between observed and solar (Fe/O) adopting (Fe/O)0 = —1.29 
from Holweger 2001). 

There are two comphcations which affect the interpretation 
of the data in Figure 7. First, if massive stars which explode 
as Type II SN also produce some Fe, the effect of time delay 
will be less pronounced for Fe than for N. Indeed, Galactic 
stars with metallicities similar to those of the DLAs consid- 
ered here, between ^ 1/10 and ~ 1/100 of solar, are defi- 
cient in Fe only by factors of 2-3. In the units of Figure 7, 
[Fe/O] = (Fe/O) - (Fe/O)© = -0.3 to -0.5 (e.g. IsraeUan et 
al. 2001). This is larger, but not by much, than the typical accu- 
racy of 0. 1-0.2 dex with which this ratio is measured in DLAs. 
A more serious problem is that, unUke O, Fe can be depleted 
onto dust and this effect will mimic the intrinsic underabun- 
dance which we are testing for. The combination of the sub- 
tlety of the effect with the uncertainty in correcting for dust de- 

* The same is not true for other Fe-peak elements, especially Zn 
whose abundance may be easier to interpret because, unlike Fe, Zn is 
not depleted onto dust. 



pletions may well explain, at least in part, why the a-element 
enhancement well established in Galactic metal-poor stars has 
proved so difficult to pin down in DLAs (Vladilo 1998; Pettini 
et al. 2000; Prochaska & Wolfe 2002, Ledoux, Bergeron, & 
Petitjean 2002). Accounting for the fraction of Fe in solid form 
requires taking into consideration the full complement of abun- 
dance measures in a DLA (e.g. Vladilo 2002), and is beyond 
the scope of this work. The values of [Fe/O] plotted in Figure 
7 should therefore be considered as indicative of the maximum 
intrinsic deficiency of Fe relative to O in the DLAs in our sam- 
ple. 

These complications make it difficult to draw firm conclu- 
sions from the results in Figure 7. Some of the DLAs appear 
to meet our expectations, by exhibiting low values of (N/O) 
matched by corresponding Fe underabundances. Similarly, 
DLAs which lie closer to the Primary line at (N/0)= —1.5 
show a range of [Fe/O] values, as expected if the time delay 
for the release of Fe is longer than that of N. But there are 
also cases which evidently do not fit this pattern. The DLA 
in Q0347— 3819 may be a special case, in that it is the only 
DLA known to date with such a strong a-element enhance- 
ment (Levshakov et al. 2002). More problematic, however, are 
the two DLAs with low (N/O) (both are upper limits) and only 
a minor (if any) Fe deficiency. Possibly we have underesti- 
mated the N abundance in these cases because of blending with 
Lyman a forest lines. ^ 

We conclude that the results of Figure 7 do not provide 
clear support for the scenario where both N and Fe lag behind 
the release of O into the ISM. Possibly, as more data become 
available, a trend between (N/O) and (Fe/O) will emerge, if one 
exists. But Figure 7 may also be telling us that there is a limit to 
the degree to which our ideas on the chemical evolution of the 
Milky Way can be applied to the high redshift DLA galaxies. 

6. Summary and Conclusions 

We have presented new measurements of the abundances of N, 
O, Si, and Fe in two metal-poor damped Lyman a systems (plus 
one sub-DLA) obtained with UVES on the VLT. We have com- 
bined these data with others from the literature to form a sam- 
ple of 10 DLAs for which reliable measurements of the abun- 
dance of N (or useful upper hmits) are available. The sample 
consists exclusively of high resolution, high S/N measurements 
obtained with 8-lOm class telescopes. Our aim was to test cur- 
rent ideas on the nucleosynthesis of N, by extending measure- 
ments of the abundances of N and O from H II regions in the 
local universe to the younger and less enriched galaxies at high 

' A similar apparent inconsistency was pointed out by Centurion et 
al. (1998) for the DLA in QOlOO+1300. The fflRES observations of 
this QSO by Lu et al. (1998) and Prochaska & Wolfe (1999), used 
here, give (N/O) + 1.5= +0.13 and [Fe/O] = -0.27, which are con- 
sistent with our working hypothesis. Centurion et al. (1998), on the 
other hand, derived a much lower colvmin density of N — by a fac- 
tor of four — from observations of the weaker N 1 AA1134 triplet with 
the 4 m William Herschel telescope on la Palma. Observations of the 
N I AA1134 triplet with an 8-10 m telescope are required to resolve 
this discrepancy. 
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redshift ({z) = 2.5) that give rise to the DLAs. Our principal 
conclusions are as follows. 

1 . The (N/O) ratio remains a difficult quantity to measure, be- 
cause the absorption lines available generally have widely 
different optical depths, so that when the N I A1200 triplet 
is detected O I A1302.2 is saturated. Furthermore, the N I 
lines are often blended within the Lyman a forest. We have 
not identified an optimum strategy for dealing with these 
difficulties. DLAs at higher redshifts than those shown here 
make accessible weaker transitions of O I at far-UV wave- 
lengths, but the chances of confusion of N I with Lyman a 
forest lines are then increased. 

2. We have argued that S, when available, can be used as a 
proxy for O, because they are both undepleted a-elements 
and the differential ionisation corrections between O I and 
S 11 should not be important when the column density 
of neutral hydrogen is as high as in DLAs. However, it 
would obviously be reassuring to confirm empirically that 
(0/S) = N(0 1)/N(S II) in at least a few cases, with spe- 
cially targeted observations. 

3. We confirm earlier reports that the (N/O) ratio shows a 
larger dispersion of values than other ratios of heavy ele- 
ments in DLAs. All 10 measurements in our sample fall 
within the region in the (N/O) vs. (O/H) plot bounded by 
the primary and secondary levels of N production. 

4. We have considered several possible explanations for this 
scatter and conclude that the most plausible one is that we 
are seeing the effects of a time delay, by about 250 Myr, in 
the release of primary N from 7 — 4 Mq stars relative to the 
massive stars which release O when they explode as Type 
II supemovae. Thus, our results provide empirical evidence 
in support of currently favoured ideas for the nucleosynthe- 
sis of N. The uniform value (N/O) ~ — L5 (on a log scale) 
seen in nearby metal-poor star-forming galaxies can be un- 
derstood in this scenario if these galaxies are not young, but 
contain older stellar populations, as indicated by a number 
of imaging studies with HST. 

5. A surprisingly high proportion (40%) of DLAs in our sam- 
ple have apparently not yet attained the full primary level 
of N enrichment at (N/O) ~ — 1.5 . Possibly, the low metal- 
licity regime — where the difference between secondary 
and primary nitrogen enrichment is most pronounced — 
preferentially selects young galaxies which have only re- 
cently condensed out of the intergalactic medium and be- 
gun forming stars. A more speculative alternative, which 
needs to be explored computationally, is that at low metal- 
licities stars with masses lower than 4M© may make a sig- 
nificant contribution to the overall N yield. The release of 
primary N may, under these circumstances, continue for 
longer than 250 Myr, perhaps for a substantial fraction of 
the Hubble time at the median (z) = 2.5 of our sample. 

6. In this scenario, the abundance of Fe — two thirds of which 
is presumed to originate from Type la supernovae — should 
also lag behind that of O, by at least as long as that of N. 
However, conflicting clues are provided by a plot (N/O) vs. 
(Fe/O) for the present data; while some DLAs do match the 
predictions of our working hypothesis, there also appear to 



be some deviant cases. The interpretation is complicated 
by the fact that some fraction of the Fe is also produced 
by massive stars and that depletion onto dust can be an is- 
sue. A larger sample of measurements is required to assess 
whether the relative abundances of N, O and Fe in DLAs 
are consistent with current ideas of chemical evolution in 
galaxies. 

7. Finally, we point out that the large dispersion in the abun- 
dance of N at these low metallicities essentially precludes 
the use of N V AA1238, 1242 doublet in constraining mod- 
els of the ionisation of the intergalactic medium from 
the ratios of different ions in the Lyman a forest (e.g. 
Boksenberg et al., in preparation). 
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Fig. 1. Normalised UVES spectra of the three QSOs observed showing the regions of the four damped Lyman a lines. Superposed on the data 
are theoretical absorption profiles for Lyman a lines with the neutral hydrogen column densities NQl 1) listed in Table 1. 
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Fig. 2. Normalised absorption profiles (histograms) of selected metal lines in the Zabs = 2.45620 DLA in Q1409+095. The j/-axis is resid- 
ual intensity. The thin continuous lines show theoretical profiles produced by VPFIT with the parameters listed in Table 2. Three velocity 
components, indicated by vertical tick marks, contribute to the absorption in this DLA. 
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Fig. 3. Normalised absorption profiles (histograms) of selected metal lines in the Zabs = 2.66820 sub-DLA in Q1409+095. The y-axis is 
residual intensity. The thin continuous lines show theoretical profiles produced by VPFIT with the parameters listed in Table 2. The fit to the 
N 1 region (within the Lyman a forest) indicates the maximum amount of N 1 which could be present. 
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Fig. 4. Normalised absorption profiles (histograms) of selected metal lines in the Zabs ~ 2.08681 DLA in Q1444+014. The y-axis is residual 
intensity. The velocity structure of this system is exceedingly complex; at least 13 separate absorption components (indicated by vertical tick 
marks in the top panel) contribute to the absorption lines. While the profile of Si II A1304.4 can be reproduced adequately by VPFIT (thin 
continuous line), we cannot derive reliable column densities for O I nor N I. In the former, the strongest components are heavily saturated, 
while the triplet structure of the latter leads to intractable blending. The VPFIT profiles for O I and N I are partial fits to the weaker components 
only, shown for illustration. 
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Fig. 5. Normalised absorption profiles (histograms) of selected metal lines in the Zabs ~ 2.07623 DLA in Q2206— 199. The y-axis is residual 
intensity. The thin continuous lines show theoretical profiles produced by VPFIT with the parameters listed in Table 2. The absorption seen in 
this DLA is due to a single velocity component with b — 6.1 km s~^. 
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Fig. 6. Abundances of N and O in extragalactic H II regions (small dots) and damped Lyman a systems (large triangles). The H 11 region 
measurements are from the studies by Kobulnicky & Skillman (1996); Thurston, Edmunds, & Henry (1996); van Zee et al. (1998); Ferguson, 
Gallagher, & Wyse (1998); Izotov & Thuan (1999); and van Zee (2000). The values for damped Lyman a systems are from the compilation in 
Table 4; filled triangles denote DLAs where the abundance of O could be measured directly, while open triangles are cases where S was used as 
a proxy for O (see text). The error bars in the bottom right-hand corner give an indication of the typical uncertainties; the large dot corresponds 
to the solar abundances of N and O from the recent reappraisal by Holweger (2001). The dashed lines are approximate representations of the 
primary and secondary levels of N production (see text). 
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Fig. 7. In this figure we plot the deficiency of N relative to the Primary level at (N/0) = —1.5 vs. [Fe/0], the departure of the (Fe/0) ratio in 
DLAs from the solar (Fe/O)© = —1.29 (Holweger 2001). In the logarithmic units used here, the abundance ratio of two elements has a typical 
measurement error of between 0.1 and 0.2 . Measured values of (N/0) and (Fe/O) for the nine DLAs shown are listed in columns (8) and (9) of 
Table 4. 



